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Abstract 

The pseudo SU(3) approach has been used to describe many low- lying 
rotational bands, as well as BE(2) and B(M1) intensities in rare earth and 
actinide nuclei, both with even-even and odd-mass numbers ||^, ^, |^, ^, ^ |^. 
The 13 P half lives of some of these parent nuclei to the ground and excited 
states of the daughter ones were evaluated for the two and zero neutrino emit- 
ting modes |Q, ^ ^, |l^, |l^ using the pseudo SU(3) scheme. The predictions 
were in good agreement with the available experimental data for ^^"Nd and 

23Sy 

The double electron capture half-lives of ^^^Dy, ^^'^Er and ^^^Yb were also 
studied using the same formalism ^]. The first nuclei was found to be the 
best candidate for experimental detection, with a half-life of « 10^* yr to the 
first excited 0+ state in ^^"^Gd. 

There are strong selection rules which restrict the two neutrino mode of 
the (3(3 decay in some other nuclei, includingjf "Gd. Experimental limits for 
the (3(3 decay of ^^"Gd have been reported Recently it was argued 

that the strong cancellation of the 2v mode in the (3(3 decay of ^'^"Gd would 
suppress the background for the detection of the Qv mode JtSt - 

In the present contribution we extend the previous research j^, ^ pH 
evaluating the 13 P half lives of ^^"Gd using the pseudo SU(3) model. While 
the 2v mode is forbidden when the most probable occupations are considered, 
states with different occupation numbers can be mixed through the pairing 
interaction. The amount of this mixing is evaluated using perturbation theory. 
The possibility of observing the /3/3 decay in ^^°Gd is discussed for both the 
2v and Qv modes. 

PACS numbers: 21.60. Cs, 21.60.Fw, 23.40. He 



1 Introduction 

The calculation of two neutrino double beta decay matrix elements has proven 
to be extremely sensitive to the details of the wave functions of the initial and 
final nuclei @. While QRPA calculations are easy to perform, the uncertain- 
ties in the residual particle-particle proton-neutron interaction strongly limit 
their predictive power . Shell model studies in the full fp shell provide reli- 
able matrix elements for **Ca and other light nuclei uM. Matrix elements 
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for Ge, Se and Xe were obtained using very large shell model spaces, 
which are however strongly truncated For heavier nuclei standard shell 
model calculations are impracticable. 

The pseudo SU(3) shell model ||^, |^ is a microscopic model which allows 
the description of heavy deformed nuclei in the laboratory frame through the 
use of a fermionic many-particle basis with good angular momentum. The 
microscopic Hamiltonian employed includes single-particle energies as well as 
pairing and quadrupole-quadrupole interactions. Most of the Hamiltonian 
parameters are fixed by known systematics. Three "rotor terms" allow a 
fine tuning of the energy spectra and are fitted for each nuclei. The pseudo 
SU(3) model has been used to describe many low- lying rotational bands, as 
well as B(E2) and B(M1) intensities in rare earth and actinide nuclei, both 
with even-even and odd-mass numbers ^, ^, ^ . It was exhibited as a 
powerful and predictive tool in the description of heavy deformed nuclei. 

The half lives of some of these parent nuclei to the ground and excited 
states of the daughter ones were evaluated for the two and zero neutrino emit- 
ting modes |Q, ^ ^, |l^ using the pseudo SU(3) scheme. The predictions 
were in good agreement with the available experimental data for ^^'^Nd and 
■^^^U. Predictions for double electron capture were also performed. These 
results are briefiy reviewed in the present contribution. 

Extending the previous research @, |, |, 0, the (3(3 half lives of ^''°Gd 
are evaluated using the pseudo SU(3) model. While the 2v mode is forbidden 
when the most probable occupations are considered, states with different oc- 
cupation numbers can be mixed through the pairing interaction. The amount 
of this mixing is evaluated using perturbation theory. The possibility of ob- 
serving the /3/3 decay in ^''"Cd is discussed for both the 2v and Qv modes. 



2 The pseudo SU(3) formalism 



In the pseudo SU (3) shell-model coupling scheme , normal parity orbitals 
(t], are identified with orbitals of a harmonic oscillator of one quanta less 
fj = rj — 1. The set of orbitals with j = j = I + s, pseudo spin s — 1/2, 
and pseudo orbital angular momentum I, define the so-called pseudo space. 
The orbitals with j — I ± 1/2 are nearly degenerate. For configurations of 
identical particles occupying a single j orbital of abnormal parity, a convenient 
characterization of states is made by means of the seniority coupling scheme. 

The many-particle states of Ua nucleons in a given shell r/a, a = v or 
TT, can be defined by the totally anti-symmetric irreducible representations 

N A 

} and {1"" } of unitary groups. The dimensions of the normal (A*') 
parity space is Q,^ = {fja -\- l){fla + 2) and that of the unique (^4) space is 
= 2{t]ci -\- 2) with the constraint ria = + ■ Proton and neutron 
states are coupled to angular momentum and in both the normal and 
unique parity sectors, respectively. The wave function of the many-particle 
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state with angular momentum J and projection M is expressed as a direct 
product of the normal and unique parity ones, as: 

I jM J2 [i-^"^ > ® iJ"" >]«■ (1) 

Since we are interested in describing low-lying energy states, only pseudo spin 
zero configurations are taken into account in the normal parity space and only 
seniority zero configurations in the abnormal parity space. This simplification 
implies that = = 0. This is a strong assumption, but one that is 
physically motivated and very useful for simplifying the calculations. 

Double beta decay, when described in the pseudo SU(3) scheme, is strongly 
dependent on the occupation numbers for protons and neutrons in the normal 
and abnormal parity states: n^,n^,n^,n^ Q. These numbers are deter- 
mined by filling the Nilsson levels from below, as discussed in Q. 

In the first series of papers Q ^, |^, |ll| we evaluated the /3/3 matrix ele- 
ments by taking into account only the leading SU(3) coupled proton-neutron 
irrep, which in recent calculation was shown to represent around 60 % of the 
wave function in even-even Dy and Er isotopes [n] 



3 Allowed two neutrino [3 /3 decays 

In all of the calculations only one active shell was allowed for protons, and 
likewise, only one for neutrons. This is a very strong truncation. For the 
Pl32v decay this implies that only one uncorrelated Gamow- Teller transi- 
tion is allowed: that which removes a neutron from a normal parity state 
with maximum angular momentum and creates a proton in the intruder shell 
(/ig/2 ^ ''•11/2 in ^^^^ earth nuclei, 111/2 ~* '''13/2 in actinides). This unique 
Gamow- Teller transition controls the /3/32v decay. Under these assumptions, if 
the occupation of the Nilsson levels is such that the number of protons in the 
abnormal states does not change for the initial and final state configurations, 
the decay is forbidden. 

The published results for the six allowed two neutrino P~ I3~ emitters are 
given in Table 1. 

Transition t'll'^lyr] T'll'^br] 

"^Nd^i'^'^Sm 2.1 X 10=*' 

148Nd ^148 g o ^ ^o20 

iso^d ^150 6.0 X 10i« 

i86w ^1*'' Os 6.1 X 10^* 5.13 X 10^5 
^^^Os ^^^^ Pt 9.0 X 10^^^ 3.28 X 10^** 

238u ^238 ^ 4 ^ -^Qll -^ 03 ^ -^q24 

Table 1: The calculated double beta half-lives for the two-neutrino and the 
zero-neutrino 

The agreement between the theoretical two neutrino half lives |Q] with the 
available data for ^^^Nd (tj/^ = 9(17) x 10^** yr) and ^^'^^U (jH'^ = 2 x 10^^ 
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yr) is good. For the theoretical /3/3oi/ half lives |8| we assumed < > = 
leV. 

4 Allowed double electron capture with two 
neutrino emission 

The 2v double electron capture (ECEC) 

2e- + (yl, Z + 2) — >{A,Z) + 2v (2) 

has larger Q- values than the concurrent 13'^ 13'^ and /?"'" EC processes, and has 
no Coulomb suppression but is very difficult to detect, because only two X 
rays are emitted together with the neutrinos. 

The double electron capture decay to excited states in the final nuclei 

{A,Z + 2) + 2e- — > {A,ZY+2v (3) 

1-^ (A,Z) + 27 

has been proposed as a good candidate to be measured The two 

gammas are far easier to detect than the X rays. A sensitivity close to ~ 
lO^^yr has been estimated for this type of experiments |p3|. 





Transition 






rV2(yr) 




^^5«Gd 


0+ ^ 


0^ 




2.74x10^^ 






0+ ^ 


0^ 




8.31x10^* 






0+ ^ 


0^ 


-(2) 


1.08 xlO^'' 


i62Er 


^le^Dy 


0+ 


0^ 




2.85x10^^ 






0+ ^ 


0^ 




3.70x10^^ 


168Yb 




0+ ^ 


0^ 




2.00x10^^ 






0+ ^ 


0^ 


-(1) 


5.36x10^^ 



Table 2: Half-lives for the ECEC2v decay to the ground and excited states 
of the final nuclei. 



In Table 2 we present the ECEC2^ decay of ^''^Dy, ^^^Er and ^'^^Yb to 
the ground and excited states of ^^"^Cd, ^''^Dy and ^^**Er respectively Q. 
The kinematical factors G2v{o') were evaluated following the prescriptions 
given in p^ . When the energy released in the decay to an excited stated 
(2_Ect) is small the available phase space G21/ is strongly reduced, and the 
half Ufe could be very large. It is the case in the double electron capture 
decay to the first excited 0^ state in '^^^Er. The nuclear matrix elements 
associated with the decay to the ground state of the final nuclei have values 
close to 0.05 - 0.06, a factor of 5 smaller than the assumption of Barabash 
|2^ , and are similar for the three nuclei studied. The nuclear matrix elements 
to excited states show a wide spread, being close to those of the ground state 
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for Dy, suppressed by a factor 5 for Er and by a factor 80 for Yb. 
While in general it is confirmed that deformed nuclei have smaller nuclear 
matrix elements than spherical j^, '^^^Dy appears to be the best candidate of 
this group for experimental detection, with a half-life around 10^* years for 
the double electron capture to the first excited 0"*" state. 



5 Forbidden two neutrino [3 (5 decays 

As mentioned above, the number of nucleons in normal and unique parity 
orbitals is determined by the filling of the deformed Nilsson orbitals. In 
this way the theory predicts the complete suppression of the I3(i2v decay for 
the following five nuclei: ^'^''Sm, "Ocd, "^^Yb, ^^^Th and ^^^'Pn [g]. It was 
expected that these forbidden decays would have, in the best case, matrix 
elements that would be no greater than 20% of the allowed ones, resulting 
in at least one order of magnitude reduction in the predicted half- life [25|. 
Experimental limits for the /3/3 decay of ^^"Gd have been reported L4|. 
Recently it was argued that the strong cancellation of the 2u mode in the /3/3 
decay of ^''"Gd would suppress the background for the detection of the Qv 
mode 

With a deformation of e = 0.26, the most probable occupations for the 
^^''Gd 14 valence protons are 8 in normal and 6 in unique parity orbitals, and 
for the 14 valence neutrons are 8 in normal and 6 in unique parity orbitals. 
The dominant component of the wave function is 

|i«°Gd,0+) = j = 0; (^13/2)^, = >a 

|{24},(10, 4)^; {24}418, 4),; 1(28, d.)K = l,J ^ >n ■ 

Assuming a slightly larger deformation for ^^''Dy, the most probable oc- 
cupations for 16 valence protons are 10 in normal and 6 in unique parity 
orbitals, and for the 12 valence neutrons are 6 in normal and 6 in unique 
parity orbitals. The dominant component of the wave function is 

|^*Dy,0+(a)> = I = 0; fe/a)'., = >a 

j{2^K(10, 4)^; {2^1418, 0).; 1(28, 4)^- = 1, J = >iv . 

The two neutrino double beta operator annihilates two neutrons and creates 
two protons with the same spatial quantum numbers. It cannot connect the 
states |^''°Gd,0+) and |^*'°Dy, 0+(a)) and the transition is forbidden. 

However, the pairing interaction allows mixing between different occupa- 
tions. In the deformed single particle Nilsson scheme it takes about Ai5 =1.50 
MeV to promote a pair of protons from the last occupied normal parity or- 
bital to the next intruder orbital. These excited state has 8 protons in normal 
and 8 in unique parity orbitals. Its wave function has the form 

|""Dy,0+(fe)) = I (/iii/2)«, = 0; (zis/a)^, = >a 

i{2*}.(10, 4),; {2='}418, 0).; 1(28, 'i)K = I, J = >n ■ 

The two neutrino double beta decay is allowed to these excited state. 
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Using perturbation theory the Dy wave function is 

r^y, 0+) = a r-Dy, 0+(a)) + 6 [""Dy, 0+(6)>, 
with + b'^ — 1 and 

b = r-Dy, 0+(fe) I Jfpai^ing I «°Dy, 0+(a)> / AE. 

Using standard SU(3) techniques we obtain a = 0.972, b — 0.233. For 
the two neutrino mode the half-Ufe is t^I^ [^^^'Dy ^^^^Gd, 0+ 0+{g.s.)] 
= l.GSxlO'^^ yr. It is delayed by one order of magnitude (compared with 
other nuclei with similar Q/j/j values) due to the presence of 6^ in the nuclear 
transition matrix element. The zero neutrino double beta transition is allowed 
to both states in ^'^"Dy due to the presence of the neutrino potential in the 

1/2 

transition operator. Assuming again < > — leV the half-life is Tq,' 
[i56Dy ^i56qj^ q+ ^ o+(ff.s.)] = 3.96xl0''5 yr @. 

The suppression of the 21/ mode would facilitate the search for the zero 
neutrino mode in the double beta decay of ^^^Dy, as it was mentioned in 
However, the value of the zero neutrino half-life is larger than the one 
expected in this reference. 



6 Conclusions 

The pseudo SU(3) shell model allows for a fully microscopic description of 
heavy deformed nuclei in the laboratory frame, in a basis having good angular 
momentum and particle number. It has proven to be a powerful and predictive 
technique to describe the low energy spectra and electromagnetic transitions 
in rare earths and actinide isotopes. It has also been applied to the description 
of the double beta decay in these nuclei. 

The model Hamiltonian includes realistic single particle energies, quadrupole- 
quadrupole and pairing interactions, whose strengths are fixed using known 
systematics. It also has three rotor terms which a allow a fine tuning of the 
energy spectra. The wave functions mix different SU(3) irreps, while in even- 
even nuclei there is always one leading irrep which contributes with about 
60% of the total wave function. 

Using only these leading irreps the double beta transition matrix elements 
were evaluated. Both the two neutrino and zero neutrino mode were studied, 
involving decays to the ground state and to excited states in the daughter 
nuclei. It was found that, at first order, the two neutrino double beta decay is 
forbidden for many candidates. The double electron capture was also studied 
for three nuclei, concluding that the ECEC of ^^''Dy could be detected in 
future experiments. 

Including difi'erent occupation numbers in normal and unique parity or- 
bitals for protons and neutrons through pairing mixing it was possible to 
evaluate the forbidden half life of ^''"Gd, which is suppressed by a factor 10 
(compared with other nuclei with similar Qpp values). It could offer a window 
to detect the zero neutrino mode which, having more decay channels, is less 
suppressed by deformation. 
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